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Abstract 

A  simple  method  is  proposed  for  battery  impedance  characterization  at  low  frequencies  (relevant  to  the  electrochemical  processes)  and  at 
higher  frequencies  (relevant  to  the  double-layer  process).  The  electrochemical  impedance  at  very  low  frequencies  is  described  using  discharge 
curve  information.  This  method  makes  it  possible  to  distinguish  between  outwardly  similar  batteries  that  differ  in  backup  times.  The  double-layer 
impedance  at  higher  frequencies  is  characterized  using  short  testing  pulses.  It  is  shown  that  the  double-layer  impedance  of  the  negative  electrode 
is  observable  in  hertz  range  while  the  positive  electrode  is  observable  at  the  boundary  to  the  millihertz  range.  The  double-layer  impedance  of 
deeply  discharged  battery  is  observable  at  higher  frequencies  than  the  impedance  of  fully  charged  (or  significantly  discharged)  battery  because  the 
electrodes  are  less  overlapping  in  frequency  and  the  spectrum  shift  to  the  higher  frequencies  is  larger  especially  for  the  negative  electrode. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  purpose  of  the  paper  is  to  propose  a  simple  method  for 
battery  impedance  characterization  at  low  frequencies  (relevant 
to  the  electrochemical  processes)  and  at  higher  frequencies  (rele¬ 
vant  to  the  double-layer  process).  The  low  frequency  impedance 
will  be  characterized  using  discharge  curves  and  the  double¬ 
layer  impedance  using  short-term  current  pulses  injected  simi¬ 
larly  to  the  direct  current  method. 

Battery  impedance  can  be  easily  measured  in  hertz  range 
[1,4].  This  provides  useful  information  on  the  double-impedance 
of  the  negative  electrode.  The  positive  electrode  should  be  anal¬ 
ysed  at  lower  frequencies  and  the  electrochemical  processes 
at  much  lower  frequencies.  The  electrochemical  impedance 
spectroscopy  [2,3]  is  designed  to  measure  slow  processes. 
These  measurements  are  time-consuming  and  limited  by  sev¬ 
eral  assumptions  like  the  time-invariance  of  a  battery  during 
the  measurements.  This  requirement  is  difficult  to  meet  at  low 
frequencies  under  discharge.  For  example,  the  measurement  at 
1  mHz  requires  17  min  of  unchanged  battery  behaviour,  which  is 
problematic  even  under  small  test  signal;  it  is  unrealistic  under 
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any  discharge.  For  example,  under  relatively  slow  discharge  of 
C//=  5  h,  the  battery  non- stationary  behaviour  is  observable  in 
the  window  of  10  s.  Also,  the  test  signal  should  be  very  small 
(to  not  corrupt  the  linear  analysis),  which  increases  the  noise 
sensitivity  of  the  measurements. 

In  this  paper,  the  proposed  method  is  not  limited  by  the  sta¬ 
tionary  or  linearity  requirements  and  is  obviously  not  granted  the 
benefits  of  linear  analysis.  The  impedance  spectrum  depends  on 
magnitude  and  form  of  the  test  signal  applied,  but  allows  for 
simple  transformation  of  a  discharge  curve  to  an  electrochemi¬ 
cal  impedance  curve  and  in  this  way  highlight  differences  of  the 
outwardly  similar  batteries  differing  in  available  backup  time. 

The  double-layer  impedances  of  positive  and  negative  elec¬ 
trodes  are  different;  they  can  be  calculated  using  another  method 
of  characterization  based  on  short-term  testing  pulses.  This 
method  allows  for  simple  transformation  of  the  injected  cur¬ 
rent  and  response  voltage  pulses  in  the  double-layer  impedance 
and  in  this  way  characterizes  the  positive  and  negative  electrodes 
separately  or  together. 

2.  Electrochemical  impedance 

Battery  electrochemical  impedance  can  be  calculated  by  the 
voltage  and  applied  current  using  a  Fourier  transformation  of 
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SOC,  % 


Fig.  1.  Battery  voltage  drop  during  discharge. 


these  processes  observed  in  certain  period.  This  idea,  justified  in 
[5],  is  simplified  in  this  paper  to  make  the  calculation  procedure 
more  practical. 

2.7.  Method 

The  discharge  curves  v[0,t]  of  two  batteries  are  shown  in 
Fig.  1  as  the  voltage  drop  v(t)  =  U(t ) — OCV  curves  of  float 
voltage  U(t)  with  respect  to  the  open  circuit  voltage  (OCV). 
They  are  relevant  to  C/Io  =  6  h  discharge  with  constant  current 
7o  =  —5  A. 

These  curves’  Fourier  transformation  is  simpler  if  the  voltage 
drop  curve  v[0,  t]  is  approximated  with  some  simpler  function, 
e.g.  an  exponential  power  series.  It  is  simpler  if  the  discharge 
curve  is  divided  into  several  pieces  and  approximated  piecewise 
as  shown  in  Fig.  2. 


SOC,  % 

100  85  70  55  40  25  10  -5 


Fig.  2.  Piecewise  approximation  of  the  voltage  drop  curve  with  model  charac¬ 
terized  with  three  sets  of  parameters. 


In  this  case,  the  second-order  exponential  power  series  is  a 
precise  enough  model 

A v(s)  =  Vo  expl^is  +  a2S2},  (1) 

where  s  is  the  time  difference,  s  =  t  —  to  with  respect  to  the 
approximation  point  to',  t  the  discharge  time,  0  <  t  <  7b;  To  the 
discharge  period  between  two  of  the  approximation  points  (pulse 
width  2h  in  Fig.  2);  A v(s)  the  voltage  drop,  A v(s)  —  v(t)  —  Vo 
with  respect  to  the  initial  voltage  Vo  =  v(to)',  obviously,  this  dif¬ 
ference  for  a  voltage  drop  v(t)  and  float  voltage  U(t)  is  the  same 
(V);  a i  and  <22  are  the  adapted  coefficients  to  the  discharge  curve. 

Theoretically,  each  piece  of  the  discharge  curve  can  be 
obtained  as  a  response  to  the  current  pulse  7o  =  5  A  inside  the 
piece  and  0  A  outside.  Thus,  the  voltage  drop  and  current  pulses 
are  absolutely  integrable  functions,  which  is  sufficient  assump¬ 
tion  for  a  Fourier  analysis. 

The  Fourier  transformations  performed  with  respect  to  the 
voltage  drop  v[to,  to  +  To]  and  current  pulse  i[to,  to  +  To]  curves 
and  transformations  divided  with  each  other  produce  simple 
result  [5] — the  following  spectral  function  for  the  electrochem¬ 
ical  impedance: 

_  f  cxp{fl|  7o  +  fl2?o  -](oTq}  _  1 

I  o\  +  2a2T0  —  )0)  a\  —  jco 

x  R]  CO 

X  1  -  exp(— jtt>7b)  ’ 

where  Zp  is  the  complex  resistance  (£2);  R  the  ohmic  resistance  at 
to,  7?=  V0//0;  7o  the  applied  current:  pulse  amplitude  (A);  To  the 
forecast  period:  pulse  width  2  h;  j  the  imaginary  unit,  j  =  T; 
co  the  angular  frequency  (rads-1),  <z>  =  27r/;/is  the  frequency 
(Hz). 

The  impedance  can  be  calculated  at  any  frequency  from  ana¬ 
lytical  formula  (2).  The  impedance  curve  is  close  to  the  real 
axis  at  high  frequencies  in  millihertz  range  because  there  are  no 
changes  in  the  electrochemical  reaction  relevant  to  this  range 
(battery  changes  slowly  during  discharge).  At  low  frequencies 
in  millihertz  range,  the  impedance  curve  is  similar  to  the  coth- 
based  Warburg  impedance  (CPE-curve  or  almost  vertical  line) 
for  positive  electrode  (or  battery),  and  to  the  tanh-based  War¬ 
burg  impedance  (semi-cycle)  for  negative  electrode.  In  general, 
the  impedance  curve  is  a  complex  exponential  function  or  log¬ 
arithmic  spiral,  spinning  in  opposite  directions  for  the  positive 
and  negative  electrodes  [5].  The  positive  electrode  is  affected 
more  severely  by  the  change  of  the  acid  concentration,  inducing 
a  counter-clockwise  spin,  while  the  spin  of  negative  electrode 
is  clockwise  oriented.  The  impedance  of  the  positive  electrode 
is  larger  than  the  impedance  of  the  negative  electrode.  Hence, 
it  dominates  in  total  for  a  battery.  It  is  most  logical  to  present 
the  impedance  curve  in  the  microhertz  range  where  most  of  the 
electrochemical  reaction  changes  take  place. 

2.2.  Results 

The  discharge  curve  in  Fig.  2  can  be  represented  as  the 
electrochemical  impedance  shown  in  Fig.  3.  It  was  calculated 
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Fig.  5.  Transient  voltage  drops  of  significantly  discharged  battery.  Measured 
and  model-predicted  voltage. 


Fig.  3.  Electrochemical  impedance  of  Battery  1  calculated  in  the  frequency 

range  1-70  |xHz  by  the  discharge  curve  in  Fig.  2.  3.1 .  Method 


by  (2)  using  adapted  coefficients  to  the  discharge  curve  of 
Battery  1. 

The  same  calculation  performed  with  respect  to  Battery  2  in 
Fig.  1  results  in  the  electrochemical  impedance  shown  in  Fig.  4. 

Battery  2  is  weaker  than  Battery  1,  which  may  also  be  con¬ 
cluded  comparing  the  discharge  curves  in  Fig.  1.  The  differ¬ 
ences  between  the  batteries  are  more  visible  by  comparing  the 
impedance  curves  in  Figs.  3  and  4  than  the  discharge  curves  in 
Fig.  1.  Therefore,  the  electrochemical  impedance  can  possibly 
distinguish  better  between  outwardly  similar  batteries  differing 
in  the  available  backup  time,  which  is  the  main  reason  to  use 
this  method  in  practice. 

3.  Double-layer  impedance 

Battery  double-layer  impedance  can  be  tested  experimentally 
using  short  current  pulses  injected  infrequently  as  in  the  direct 
current  method. 


Fig.  4.  Electrochemical  impedance  of  a  weaker  Battery  2. 


The  transient  voltage  induced  by  the  imposed  square-wave 
current  pulses  can  be  represented  with  a  two-rate  differential 
model  for  the  positive  and  negative  electrodes  [2].  This  model 
solved  analytically  [2]  can  be  represented  as  the  sum  of  two 
exponential  functions 

j£)  -  * exp 

(3) 

where  s  is  the  time  difference,  s  =  t  —  to  with  respect  to  the 
approximation  point  to',  t  the  testing  time,  0<t<T\;  T\  the 
testing  period:  pulse  width  2  s  in  Figs.  5  and  6;  A u(s)  the  volt¬ 
age  drop,  Au(s)  =  u(t)  —  v(to)  with  respect  to  the  initial  drop 
v(to)  =  IoR  (approximation  point)  (V),  the  voltage  drop  u(t)  is 
defined  as  a  response  of  voltage  to  the  applied  dc  current  Io  and 
to  the  injected  ac  current  pulse  A/;  /  the  overall  applied  current, 
I = Io  +  AI  (A);  R  the  overall  ohmic  resistance,  R  =  Ro  +  R\  +  R2, 
R\ — charge-transfer  resistance  of  positive  or  R2 — negative  elec¬ 
trodes,  Ro — total  resistance  of  conducting  elements,  it  includes 
resistance  of  electrodes  and  separator,  grid  resistance,  intercell 
connector  resistance  and  battery  terminals  resistance  (£2);  C\  is 


A  u(s)  =  AI\ 


R  —  R\  exp 


Fig.  6.  Transient  voltage  drops  of  deeply  discharged  battery.  Measured  and 
model-predicted  voltage. 
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Table  1 


Model  parameters:  estimated  by  measurements 


SOC  (%) 

Ro  (mf2) 

Positive  electrode 

Negative  electrode 

R\  (m£2) 

Ci  (F) 

R2  (mfi)  C2  (F) 

9 

2 

19 

8 

11  0.2 

27 

2 

13 

16 

7  8 

the  double-layer  capacitance  of  positive  or  C2 — negative  elec¬ 
trodes  (F). 

These  parameters  can  be  calculated  as  mean  squares  esti¬ 
mates,  fitting  the  model  (3)  with  a  measured  voltage  drop  curve. 
They  depend  on  the  state  of  charge  (SOC)  and  are  different 
for  significantly  discharged  (SOC  =  27%)  and  deeply  discharged 
(SOC  =  9%)  battery  (Table  1). 

Their  model  fit  is  shown  in  Figs.  5  and  6,  correspondingly. 

The  double-layer  impedance  can  be  calculated  as  a  Fourier 
transformation  Zdl  =  F{Am(^)  x  AI~[(s)}  of  the  voltage  drop 
curve  and  current  pulse  curve.  The  resulting  impedance  can  be 
represented  as  the  following  spectral  function: 


Zdl  =  R  + 


R 1 

1  T-  j  coR\  C\ 


+ 


Ri 

1  -\-iojR2C2 


The  inductive  components  were  not  included  in  model  (3)  and 
not  represented  in  spectrum  (4).  For  this  reason,  the  method  is 
not  applicable  at  higher  frequency  than  1  kHz.  There  may  be  also 
some  inaccuracy  in  impedance  spectrum  at  frequencies  close  to 
this  limit. 


3.2.  Results 

The  double-layer  impedance  of  the  significantly  discharged 
(SOC  =  27%)  battery  is  shown  in  Fig.  7;  it  was  calculated 
by  (4)  for  the  positive  and  negative  electrodes  and  in  total 
for  the  battery.  Impact  of  the  positive  electrode  on  the  total 
battery  impedance  is  larger  than  the  impact  of  the  negative 
electrode. 


Fig.  8.  Modulus  of  significantly  discharged  battery,  shown  for  the  positive  and 
negative  electrodes  and  in  total. 


Fig.  9.  Phase  angle  of  significantly  discharged  battery,  shown  for  the  positive 
and  negative  electrodes  and  in  total. 


Fig.  7.  Double-layer  impedance  of  significantly  discharged  battery,  shown  for  Fig-  10-  Double-layer  impedance  of  deeply  discharged  battery  shown  for  the 
the  positive  and  negative  electrodes  and  in  total  for  the  battery.  positive  and  negative  electrodes  and  in  total. 
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Fig.  11.  Modulus  of  significantly  discharged  and  deeply  discharged  batteries. 


Fig.  12.  Phase  angles  of  significantly  discharged  and  deeply  discharged  batter¬ 
ies. 

This  double-layer  impedance  is  also  represented  as  a  spec¬ 
trum  of  modulus  and  phase  angle  in  Figs.  8  and  9. 

The  double-layer  impedance  of  the  negative  electrode  is 
clearly  observable  in  hertz  range;  for  the  positive  electrode  and 
total  impedance,  it  is  observable  at  lower  frequencies  as  shown 
in  Figs.  8  and  9;  for  fully  charged  battery,  it  is  observable  at  the 
boundary  to  the  millihertz  range  or  lower. 

The  behaviour  of  the  deeply  discharged  (SOC  =  9%)  battery 
is  different  from  the  significantly  discharged  and  fully  charged 
batteries.  The  double-layer  impedance  is  larger  and  the  positive 


and  negative  electrodes  are  less  overlapping  in  frequency,  which 
produces  the  double  arc  in  the  overall  battery  curve  in  Fig.  10. 

There  is  also  a  large  shift  of  the  modulus  and  phase  angle 
to  the  higher  frequencies  for  the  deeply  discharged  battery 
because  spectrum  shift  is  larger  for  the  negative  electrode  in 
Figs.  11  and  12. 

4.  Battery  overall  impedance 

Battery  overall  impedance  of  the  electrochemical  and  double¬ 
layer  processes  can  be  calculated  approximately  by  the  super¬ 
position  of  both  impedances 

Z  =  Zp  +  Zdl,  (5) 

where  Z  is  the  overall  impedance;  Zdl  the  double-layer 
impedance;  Zp  is  the  electrochemical  impedance. 

This  approximation  is  justified  by  the  nature  of  slow  and  fast 
processes.  They  are  not  overlapping  in  frequency  [5] :  the  double¬ 
layer  capacitance  is  an  effective  shunt  for  the  Faradic  current  at 
the  higher  frequencies  than  0. 1  Hz. 

5.  Conclusion 

It  was  shown  that  a  battery  can  be  characterized  at  very  low 
frequencies  with  the  electrochemical  impedance  calculated  from 
the  discharge  curve  and  at  higher  frequencies  with  the  double¬ 
layer  impedance  calculated  from  the  voltage  response  pulses  to 
short  current  pulses.  The  impedance  characterization  may  give 
better  information  on  the  battery’s  internal  capability  to  perform 
than  the  discharge  curve  directly.  It  may  also  be  possible  to  char¬ 
acterize  the  battery  performance  using  shorter  discharge  period 
(SOC  >  70%)  than  the  complete  capacity  test  when  model-based 
methods  are  utilized. 
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